Pions and kaons from stringy quark matter 



T. S. Biro and K. Urmossy 

MTA KFKI Research Institute for Particle and Nuclear Physics 
H-1525 Budapest, P.O.Box 49. 

E-mail: tsbiroOrmki.kfki.hu, ukarolyOfreemail.hu 
PACS numbers: 21.65.Qr, 25.75.Ag 



Abstract. Different hadron transverse momentum spectra are calculated in 
a non-extensive statistical, quark-coalescence model. We suggest that extreme 
relativistic kinematics may be responsible for power-law tails under quite general 
assumptions and a quark scaling could be recognized in the powers. For the low 
energy part also a gluonic string contribution is needed to describe experimental 
data: its length distribution and fractality are fitted. 



1. Introduction 

The hadronization of quark matter in relativistic heavy ion collisions is a key 
phenomenon for understanding the interplay between quark confinement and 
deconfincmcnt, between statistical and dynamical factors. The purely hadronic 
statistical model [TJ [3J [H [5] was successful in describing several number ratios of 
different hadron species and - by assuming a collective flow - some bulk properties of 
individual hadron spectra. This model fits experimental results at several bombarding 
energies ranging from SIS to RHIC, and finds that the main parameters, the 
temperature T and the baryochemical potential, /i lies in a narrow range near to 
a line described by a constant energy per particle of E/N = 1 GeV. This value is 
six times larger than the apparent temperature of spectra after subtracting transverse 
flow effects. This raises the question that how can be a quark matter, describable 
by perturbative QCD and containing nearly massless quarks and massless gluons, 
transformed into a hadronic matter with an average mass of m = 750 MeV (based on 
the simple nonrelativistic ideal gas formula E/N = m + 3T/2). 

Recently we have suggested an approach to massless quark matter supplemented 
with stringy interaction corrections [7] giving a contribution to the free energy density 
which is proportional to the (colored) number density on a fractional power. It was 
shown that such an equation of state accommodates the E/N = 6T ratio at zero 
pressure, at the edge of mechanical stability. This model assumes underlying non- 
perturbative effects around and above the color deconfinement temperature, T c « 170 
MeV. Quark matter equation of state data from lattice QCD support an interaction 
modification to the ideal massless gas picture of /i„t oc n 2 / 3 . Since the basic concept 
is quite simple, we repeat briefly this argumentation in the present paper. 

On the other hand the quark coalescence picture proved to be applicable to 
the hadronization problem in relativistic heavy ion collisions at RHIC energies in 
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several respects. The quark number scaling of asymmetric azimuthal elliptic flow 
factor [T51 [T^l [T31 IT51 1 1 7j . v%, and the number ratios of heavy flavor, most prominently 
strange, hadrons to the light ones, as it has been repeatedly reported in calculations 
using the ALCOR model [18], all underline the viability of this concept. We also wish to 
add a further observable, namely the Tsallis parameter, q— 1, regulating a cut power- 
law fit to transverse momentum spectra, to this list. In this paper we demonstrate 
that this parameter also may follow a quark number scaling (coalescence) law. 

In order to interpret this result in proper context, also some introduction 
is necessary into recent developments on the interpretation of non-extensive 
thermodynamics. We give a comprised summary of investigating the thermodynamical 
limit of abstract composition rules[8j for two particle systems, which motivates the 
use of cut power-law distributions of quark momenta in the framework of statistical 
physics models. By doing so we connect predictions from dynamical parton models 
for pT-spectra to those of simpler statistical models smoothly |12j. 



2. String distribution and equation of state in quark matter 

Assuming a pair interaction energy also in the quark matter, which is proportional 
to the nearest neighbour distance, one deals with an energy density modification of 
no (£) = An 2 /*. The free energy density becomes 

f(n,T) = f id (n,T)+An 2 / 3 . (1) 

Deriving the energy density and pressure via standard thermodynamical formulas one 
obtains density dependent, stringy mean field terms, which - unlike in the bag model 
- do not compensate each other exactly: 

e = e id + \An 2 ^, p = p id - ^An 2 ' 3 . (2) 

At vanishing pressure the correction can be related to the ideal pressure by p = and 
the energy per particle for massless constituents in a Boltzmann gas (pid = nT) is 
given by 



E_e 

N ~ n 



e id + 3p i( 
p=o n 



^ = 6T. (3) 



We note that any interaction correction proportional to n 1 7 in the free energy leads 
to the following scaling of the interaction measure at high temperature: 

^ ~ T- 3 -- 1 (4) 
Lattice QCD calculation results definitely favor 7 = 1/3. 



3. Power-law tailed distribution from repeated composition rules 

Power-law tailed distributions are numerous in physical phenomena, they can be 
viewed as stationary distributions under certain constraints. A general approach 
in the framework of non-extensive thermodynamics has been presented in Ref.JSJ- 
Starting with a general non-additive composition rule, its manifold repetition becomes 
asymptotically an associative rule. Among these statistical attractors the rule 
h(x, y) = x + y + axy is particularly simple (here the simple addition is reconstructed 
for a = 0). It is asymptotic to any rule of h(x,y) — x + y + G(xy) type. We have 
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shown that extreme relativistic kinematics can lead to such a rule for kinetic energy 
composition if the pair interaction is a function of the relativistic momentum transfer 
squared, Q 2 , only. We surmise that parton cascades based on pQCD cross sections in 
the large collision number limit may belong to this class of effective composition rules. 

Any associative composition rule can be expressed by the addition of formal 
logarithms. In our case (1 + ah) = (1 + ax)(l + ay) so the formal logarithm is given by 
L a (x) = (1/a) ln(l+aa;). The canonical energy distribution by this energy composition 
rule is given by 

f cq (E) = e-*** = (1 + aE)^ aT . (5) 

This power-law tailed energy distribution is the Pareto-Tsallis distribution with 
q = 1 + aT. Its reciprocal logarithmic slope is proportional to the energy: 

T(E)=T+(q-l)E. (6) 



4. Coalescence quark scaling in the power-law 

Besides the familiar factorizing coalescence formula valid at high energies and 
momenta, we consider a hadronization kinematics with string like objects contributing 
to the total energy of hadrons. Regarding an energy balance of 



E t + E 2 + a£ = yjm% + (pi + p 2 ) 2 , (7) 
for meson formation (and a three-fold one for baryon formation) we consider 

f(E H ) = Jd 3 pdE 1 d 3 qdE 2 dmf 1 + ?j f 2 Qp - ?j 

C{m, q) g(m) 9 (E 2 H - (E 1 + E 2 + mf) (8) 

for the coalesced hadron spectra. Here Q(x) denotes the step function, and the 
integration variable m = a£ runs over the energy values put into a string. C(q,m) 
is the coalescence factor, about which here we assume that prefers very low relative 
momenta only and put C(q,m) = cS 3 (q). About the string length distribution we 
assume the form 



g(m) r*> m exp 



r(i + i/d)-HL 

(m) 



(9) 



argued for in Ref.[19 . Here the effective dimensionality of the moving string, d, and 
the average energy stored in it, (m) are parameters. The coalescence ratio can then 
be interpreted in terms of the string mass distribution: 

f H (E H ) 



fi(Ei)f 2 (E 2 



g{m)dm. (10) 

Eh-E!-E 2 



5. Results 

In the following figures thermal Tsallis-Pareto distributions are fitted to several 
hadronic transverse momentum spectra obtained at RHIC by the PHENIXpOl ETj 
and STAR[22j |23l [M [25] groups. In Fig[T] the differential yields for pions, kaons, 
eta and phi mesons, protons, Lambda-s and Xi hyperons are shown including their 
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Figure 1. Tsallis distribution fit to different hadronic transverse momentum 
spectra with a transverse blast wave taken into account. 
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Figure 2. Tsallis distribution fit parameters for quark matter extracted from 
hadron spectra. The slope scales linearly with the mass due the Tsallis formula 
and follows the valence quark number scaling indicated by the red and blue 
lines. The black line belongs to the conjectured Tsallis-Pareto distributed quark 
matter, the star denotes the average hadron mass value of 750 MeV extracted 
from E/N = m + 3T/2 = 1 GeV at T c = 175 MeV temperature. 



antiparticle partners. The red continuous lines represent Tsallis-Pareto fits with best 
parameters for each hadron. Transforming these spectra in a frame co-moving with 
a transverse blast wave, the co-moving energy distribution f(E) is given in terms of 
E = 7(771^ — vpx)- Since this transformation is mass dependent, the pt spectra show 
maxima at different values p™ ax = jvm. 

We assume that already the quark matter had this collective transverse flow at 
the hadronization, only the hadron spectra were formed according to the coalescence 
formula eq.(fT0]). FigO plots the extracted quark matter q- value from different hadron 
spectra (upper part). This value is quite stable around q m 1.2. The slope parameters 
(lower part) at the minimal co-moving energy, E = m, (meaning px = "fvm, the 
position of the maxima) are proportional to the hadron mass due to eq.©. This 
proportionality can be observed even after disentangling the transverse flow effect, in 
contrast to non-Tsallis (q = 1) fits done earlier to heavy ion data. 

The factorization of hadron spectra to the respective quark spectra, as predicted 
by the naive coalescence model, is well satisfied at high momenta. For lower momenta, 
however, this scaling is violated. It can be seen by the ratio /jj/Zi/a shown in FigJ3] 
(left part). 

In order to explain this violation we supplement the quark coalescence picture by 
that of strings built in into the forming hadrons. Assuming that all strings longer 
than a threshold value contribute with uniform probability to the corresponding 
hadron yield, the missing energy in the coalescence model can be attributed to 
strings, distributed in their length. Fig[3] (right part) plots the differential string 
energy distribution g(m). Based on experimental pion, kaon and antiproton spectra 
in RHIC STAR and PHENIX experiments, we fit the average mass contribution values 
( m ) OT = 60 MeV, (m) K = 120 MeV, (m) F = 185 MeV and dimensionalities = 1.30, 
d K = 1.60 and dp = 1.67. 
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Figure 3. Violation of the coalescence product rule at low energy (left) Prediction 
for the differential string-energy distributions by coalescence from Tsallis-Pareto 
distributed quark matter and strings (right). 
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